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ABSTRACT 

RXJ1856. 5-3754 and RX J0720.4-3125 are the only young isolated radio-quiet neutron stars 
(NSs) for which trigonometric parallaxes were measured. Due to detection of their thermal 
emission in X-rays they are important to study NS cooling and to probe theoretical cooling 
models. Hence, a precise determination of their age is essential. 

Recently, new parallax measurements of RXJ1856. 5-3754 and RXJ0720.4-3125 were ob- 
tained. Considering that NSs may originate from binary systems that got disrupted due to an 
asymmetric supernova, we attempt to identify runaway stars which may have been former 
companions to the NS progenitors. Such an identification would strongly support a particular 
birth scenario with time and place. 

We trace back each NS, runaway star and the centres of possible birth associations (assuming 
that most NSs are ejected directly from their parent association) to find close encounters. The 
kinematic age is then given by the time since the encounter. We use Monte Carlo simulations 
to account for observational uncertainties and evaluate the outcome statistically. 
Using the most recent parallax measurement of 8. 16 ±0.80 mas for RX J 1856.5-3754 by Wal- 
ter et al., we find that it originated in the Upper Scorpius association 0.46 ± 0.05 Myr ago. 
This kinematic age is slightly larger than the value we reported earlier (0.3 Myr) using the old 
parallax value of 5.6±0.6 mas by Kaplan. Our result is strongly supported by its current radial 
velocity that we predict to be 6I20 km/s. This implies an inclination angle to the line-of sight 
of 88 ± 6 deg consistent with estimates by van Kerkwijk & Kulkarni from the bow shock. No 
suitable runaway star was found to be a potential former companion of RX J 1856.5-3754. 
Making use of a recent parallax measurement for RXJ0720.4-3125 of 3.6 ± 1.6 mas by 
Eisenbeiss, we find that this NS was possibly born in Trumpler 10 0.85 ± 0.15 Myr ago. 
This kinematic age is somewhat larger than the one obtained using the old parallax value 
of 2.77 ± 1.29 mas by Kaplan et al. (0.5 Myr). We suggest the BO runaway supergiant HIP 
43158 as a candidate for a former companion of the progenitor star. Then, the current distance 
of RX J0720.4-3125 to the Sun should be 286^23 P*^' ^'^ agreement with recent measurements. 
We then expect the radial velocity of RX J0720.4-3125 to be -76+^7 km/s. 
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1 INTRODUCTION 

Neutron stars (NSs) show large proper motions which, with known 
distances, indicate high space velocities (e.g. Lyne & Lorimer 
1994; Lorimer et al. 1997; Hansen & Phinney 1997; Cordes & 
Chernoff 1998; Arzoumanian et al. 2002; Hobbs et al. 2005). 
Some NSs even show velocities of the order of ~ 1000 km/s (e.g. 
PSR B1508+55, Chatterjee et al. 2005; PSR B2223-^65, Harrison 
et al. 1993; Taylor & Cordes 1993; RX J0822-4300, Hui & Becker 
2006, Winkler & Petre 2007). Those high velocities may be the 
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result of asymmetric supernova (SN) explosions assigning the 
new-born NS a kick velocity for that a number of mechanisms have 
been suggested (e.g. Burrows & Hayes 1996; Janka & Mueller 
1996; Janka et al. 2005; Wang et al. 2006; Kisslinger et al. 2009). 
Another possibility is that the high-velocity NSs are the remnants 
of (symmetric^) SN explosions of the so-called hyper- velocity 



^ Although numerous three-dimensional simulations showed that SNe are 
most likely asymmetric (e.g. Wongwathanarat et al. 2010; Hammer et al. 



2010; Dimmelmeier et al. 2008; 
Blondin & Mezzacappa 2007). 



Janka et al. 2008; Scheck et al. 2006; 
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runaway stars (Gvaramadze 2007; Gvaramadze et al. 2008; 
Gvaramadze 2009) which were ejected due to dynamical three- or 
four-body encounters either from the Galactic Centre (Hills 1988) 
or from massive star clusters in the Galactic disk. 
About 30 % of young stars show different velocity properties than 
normal Population I stars (Stone 1991; Blaauw 1993; Tetzlaff 
et al. 2011, thereafter Til, a catalogue of young runaway stars). 
Two scenarios are accepted to produce those so-called "runaway 
stars" (Blaauw 1961). The binary-SN scenario (Blaauw 1961) is 
related to the formation of the high velocity NSs (but note that 
high-velocity NSs may also be the result of an asymmetric SN 
explosion of a single massive star or a SN of a massive runaway 
star, see above): The runaway and NS are the products of a SN 
within a binary system. The velocity of the former secondary 
may be as large as its original orbital velocity (Tauris & Takens 
1998). Runaway stars produced in this scenario should share 
typical properties such as a high rotational velocity vsini and 
an enhanced helium abundance owing to momentum and mass 
transfer during binary evolution (Blaauw 1993). We refer to such 
runaway stars as BSS runaway stars. The second scenario is 
the dynamical ejection due to gravitational interactions between 
massive stars in dense clusters (Poveda et al. 1967). 
Several studies have been made to investigate the origin of runaway 
stars (e.g. Blaauw 1961; Gies & Bolton 1986; Hoogerwerf et al. 
2001; de Wit et al. 2005; Schilbach & Roser 2008; Gvaramadze & 
Bomans 2008) but only a few regarding the origin of fast-moving 
NSs (e.g. Hoogerwerf et al. 2001; Vlemmings et al. 2004; Tetzlaff 
etal. 2010, thereafter TIO). 

In this paper we will re-investigate the origin of two isolated 
radio-quiet X-ray emitting NSs, namely RX J 1856.5-3754 and 
RX J0720.4-3125. So far, only seven such sources have been 
identified for which they were named "The Magnificent Seven" 
(M7) (Treves et al. 2001; for recent reviews see Haberl 2007; 
Kaplan 2008). They are bright X-ray sources associated with 
faint blue optical counterparts arising from a hot cooling surface. 
In only the two cases studied in this paper a trigonometric 
parallax was obtained. The M7, and especially RXJl 856.5-3754 
and RXJ0720.4-3125, are very important NSs since brightness, 
parallax and temperature yield the size of the emitting area and 
hence, their (model dependent) radii. From spectra, one can in 
principle determine their masses and atmosphere composition, 
which eventually may lead to constraints on the equation of state. 
With known luminosity and age, cooling curves can be verified. 
As the characteristic age represents only a rough estimate of the 
true age (e.g. Blandford & Romani 1988; Gaensler & Frail 2000; 
Migliazzo et al. 2002), the kinematic age is very important to 
have a better estimate of the true age. Also, the characteristic age 
is significantly influenced by pulsar winds (e.g. Wu et al. 2003) 
and possibly by emission of gravitational waves (e.g. Wette et al. 
2008). 

Under the assumption that most NSs form in associations or 
clusters of massive stars, the identification of the birth association 
of a NS by its flight path is possible if the uncertainties in the 
distance are moderate. Our assumption is justified since we observe 
associations and clusters of massive stars (also, the dispersion time 
scale of a massive star cluster is much longer than the lifetime of 
its most massive member stars), i.e. only a small fraction of their 
massive member stars is ejected from their parent cluster due to 
gravitational interactions before they end their lives in SNe (about 
20-30 per cent of O and early B type stars are located outside of 
clusters. Mason et al. e.g. 1998; Maiz-Apellaniz et al. e.g. 2004, 
i.e. >70 per cent of a cluster's member stars remain in the cluster). 



There are certainly NSs that form outside their parent cluster as 
there are massive runaway stars that will later explode in a SN 
event and become NSs. For the latter, the identification of their 
parent associations or birth sites is hardly possible. Also, it is 
possible, that a massive binary was ejected from its parent cluster. 
If the primary experiences a SN, the secondary might be ejected 
from the system and will later explode in a SN. For NSs whose 
progenitors experienced this so-called two-step-ejection scenario 
(Pflamm-Altenburg & Kroupa 2010), it is not possible to identify 
their formation sites. However, this scenario applies only to 1-4 
per cent of O stars (Pflamm-Altenburg & Kroupa 2010). NSs 
formed from runaway stars may show a higher space velocity than 
those formed in their parent cluster since their velocity vector is a 
superposition of the runaway star's velocity and the kick velocity 
the NS receives at birth. 

When searching for the parent association of a NS (under the 
above assumptions), the result is often not unique (TIO). For that 
reason, it is desirable to find a second indicator for the location and 
time of the past SN, e.g. a possible former companion star that is 
now a runaway star 

Here, we note that not for every NS a (identifiable) runaway 
star must exist. The NS progenitor could have been a single 
star (possibly a runaway star, see above); a former massive 
companion could have already undergone a SN and is now a 
NS; or the former companion that is now a runaway star was not 
yet identified as such because, e.g., its absolute velocity is low 
(e.g. Tauris & Takens 1998) or the direction of its velocity vec- 
tor is not significantly different from those of its neighbouring stars. 

After investigating possible parent associations for 
RXJ1856.5-3754 and RX J0720.4-3125 in sections 3.1 and 
3.3, respectively, applying the procedure described in section 2 
and using most recent parallax measurements, we attempt to 
identify the possible former companion for RXJ1856. 5-3754 
and RXJ0720.4-3125 in sections 3.2 and 3.4, respectively. We 
summarise our results and draw our conclusions in section 4. 



2 PROCEDURE 

Basically, we applied the same procedure as already carried out 
in TIO, to that we refer for details (details may also be found 
in the additional online supporting information'^). We perform 
Monte Carlo simulations to account for the uncertainties of the 
observables and the unknown radial velocity of NSs for which we 
assume a probability distribution derived from the one of pulsar 
space velocities by Hobbs et al. (2005). 

In this paper, we extended the sample of 140 young as- 
sociations and clusters given in TIO to a sample of 295 young 
associations and clusters within 3 kpc from the Sun. This sample 
of young associations and clusters contains 140 associations and 
clusters listed in TIO, four additional young nearby associations 
from Torres et al. (2008) (Columba, Carina, Octans, Argus), 101 
clusters with at least one star with a spectral type earlier B3 in 
the WEBDA database'* as well as 50 additional young clusters 
with kinematic data available in cluster catalogues (WEBDA; 



^ http://www.astro.uni-jena.de/~nina/supporting_info.pdf 

■* http://www.univie.ac.at/webda/webda.html; Mermilliod & Paunzen 

(2003). 
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Kharchenko et al. 2005, 2007; Dias et al. 2010). 

The distribution of separations dmin (that is the smallest 
separation between the NS and the association/ cluster centre 
found in each Monte Carlo run) is supposed to obey a distribution 
of absolute differences of two 3D Gaussians (equation 1 in TIO; 
equation A3 in Hoogerwerf et al. 2001; see also additional online 
supporting material). For investigating encounters with runaway 
stars, we calculate dmin between the NS and the runaway star. If 
the two stars were at the same time and place, the dmin distribution 
will follow equation 2 in TIO (equation A4 in Hoogerwerf et al. 
2001). 

We will adapt the theoretical formulae only to the first part of the 
dmin distribution such that the slope and peak can be explained. 
The parameter /i then gives the positional difference between the 
two objects. The error of this difference. A, can be estimated as 
A^ = af + ai = 2(7^ 

To associate the NS/runaway star encounter position with 
an association/cluster, the trajectory of the latter is calculated 
simultaneously. Runaway star data are taken from Til (and 
references therein, mainly van Leeuwen 2007) for 2547 runaway 
stars (1705'* with full 3D kinematics). 

All calculations of trajectories are performed in a co- 
ordinate system centred on the Sun at present. We ac- 
count for Solar motion using a local standard of rest of 
{UVW)q = (10.4 ±0.4, 11.6 ±0.2, 6.1 ±0.2) km/s (TU), 
where U, V, W are the velocity components in right handed 
Cartesian coordinates. 

In general, we first perform 10** Monte Carlo runs for 
each NS/association pair and 10'* Monte Carlo runs for each 
NS/runaway star pair (for the latter less runs are still sufficient 
due to the smaller errors on the runaway star kinematics compared 
to the dispersion of the association velocities) to find those asso- 
ciations and runaway stars that potentially crossed the past path 
of the NS, i.e. those for which the smallest dmin value found 
in the calculations is less than three times the association radius 
(for NS/association pairs) or less than 10 pc (for NS/runaway star 
pairs), respectively. The latter value is justified since even smaller 
dmin values are expected to be found after lO"* runs if the NS and 
the runaway star once were at the same place (see additional online 
supporting information). Those associations and runaway stars that 
fulfilled those conditions are then selected for a more detailed in- 
vestigation (one to three million Monte Carlo runs). The outcome 
of these simulations is then discussed in detail. Regarding associa- 
tions, we search for those for which the NS could have been within 
the association boundaries in the past while for runaway stars, we 
are looking for those runaway stars for which the NS and the run- 
away star might have been at the same place in the past, hence the 
distribution (slope) should obey equation 2 in TIO. After three mil- 
lion runs the smallest dmin value found is expected to be smaller 
than one parsec (see additional online supporting information). If 
this criterion is satisfied, we adapt the theoretically expected dis- 
tributions to the first part of the "observed" dmin distribution to 



explain its slope and peak. If we find /i ~ 0, we select the runaway 
star to be a former companion candidate. 



3 RESULTS 

3.1 Identifying the parent association of RX J1856.5-3754 

Walter et al. (2010) recently reported a parallactic distance of 
RXJ1856. 5-3754 of 123l}5pc confirming earlier measurements 
by Walter & Lattimer (2002). This distance is significantly smaller 
than the value of 178 ± 20 pc claimed by Kaplan (2003) and used 
by TIO to evaluate the birth place of RXI1856.5-3754. Hence, it is 
worthwhile to re-investigate the origin of RX 1 1856.5-3754. 
We adopt the following parameters for RXJ1856. 5-3754 for the 
right ascension a, declination 5, parallax n and proper motion p 
(Walter etal. 2010): 



Q = 18''56'"35?795, 5 = -37°54'35"54, 

TT = 8.16 ±0.80 mas, 

Ha = 325.9 ± 2.3 mas/yr, 

/if = —59.2 ± 2.1 mas/yr, 



(1) 



where fi'^ is the proper motion in right ascension corrected for 
declination. 

Given its proper motion and parallax, the transverse velocity 
of RXJ1856.5-3754 is vt = I92I2I km/s. van Kerkwijk & 
Kulkarni (2001) tried to measure the inclination i of the bow shock 
which RXJl 856.5-3754 creates in the interstellar medium (ISM). 
Unfortunately, they could not obtain a precise result but estimated 
i = 60 ± 15 deg or i even closer to 90 deg, depending on the 
model. A lower limit on i of 45 dog implies a maximum radial 
velocity modulus of ~ 250 km/s (3a Vt). We will address this 
issue later. 

First, we perform lO'* Monte Carlo runs to find close en- 
counters between RX 1 1856.5-3754 and any association/cluster in 
the past five million years'. We select those associations/clusters 
for which the smallest separation dmin found was less than three 
times the association/cluster radius, 18 in total. For those 18 as- 
sociations/clusters, we carry out another one million Monte Carlo 
runs. For 13 of them we find close encounters consistent with 
the association/cluster boundaries; however, most of them [Tu- 
cana/Horologium (Tuc-Hor), the /3 Pictoris group (/3 Pic-Cap), 
AB Doradus (AB Dor), Hercules-Lyrae, Sagittarius 0B5, Scor- 
pius 0B4 (Sco 0B4), Pismis 24, Trumpler 27, NGC 6383, van den 
Bergh-Hagen 217, NGC 6396] can be excluded because in these 
cases the radial velocity would need to be \vr\ > 350 km/s (larger 
than maximum Vr inferred from bow shock, cf. above discussion 
on the bow shock). 

We adapt the theoretically expected distribution to the first part of 
the dmin distribution such that the slope and peak can be adjusted 
(see also Fig. 1 for an example). Also, in all but one (Sco 0B4) 
of these cases, adapting the theoretical curve to the dmin distribu- 
tion suggests that the closest approach of RX J1856.5-3754 to the 
association/cluster was outside the association/cluster boundaries. 
For the extended Corona-Australis association (Ext. R CrA), small 
separations between RXJ1856. 5-3754 and the association centre 



■* In Til, 1703 stars are listed with radial velocity (Vr) measurements. For 
two further stars, we find Vr in the literature: HIP 17158 with Vr = 55.19± 
0.73 km/s (Javakhishvili & Salukvadze 1995) and HIP 68281 with Vr = 
-11 km/s (Stetson 1983). 



"" Comparing their effective temperatures with cooling curves, RX Jl 856.5- 
3754 and RX J0720.4-3125 cannot be older (see e.g. Page et al. 2004; Lat- 
timer & Prakash 2004). 
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are found for 0.04 ± 0.02 Myr in the past. This is not surprising 
because the present position of RXJ1856. 5-3754 lies within the 
association and it would need up to ~ 0.2 Myr to cross the as- 
sociation (assuming a maximum space velocity of 350 km/s, i.e. 
3cr Vt and Vr,max = 250 km/s, and an extension of Ext. R CrA of 
62 pc). 

Hence, Ext. R CrA is not the parent association of RXJ1856.5- 
3754 unless the NS is only a few 10000 yr old, however then we 
would probably still see the SN remnant (SNR) (there is no SNR 
known in this area. Green 2009; A. Poghosyan, priv. comm.). The 
only association that fits perfectly to being the parent association 
of RX J 1856.5-3754 is Upper Scorpius (US) (as suggested before 
by Walter & Lattimer 2002 using a much simpler calculation and 
TIO using a larger distance). Table 1 (a) summarises the derived 
NS parameters and the position and time of the SN (for the de- 
duction of the properties see TIO or additional online supporting 
information). The displacement between RXJ1856. 5-3754 and the 
US centre would be 7.9 pc as inferred from the theoretically ex- 
pected curve. The dmin distribution and the distribution of the cor- 
responding flight times r are shown in Fig. 1, top panel. The SN 
would have occurred well within US that has a radius of ~ 15 pc. 
We find that the radial velocity of RXJl 856.5-3754 is probably 
very small (now indicated by both the parent association and the 
bow shock) although utilising the velocity distribution for pulsar 
velocities gives priority to larger values. 

Since both the radial velocity and the parallax are strongly cor- 
related as a larger distance can be compensated for by a larger 
radial velocity, parallaxes may be slightly biased towards smaller 
values (since larger radial velocities are preferred in the calcula- 
tions). Therefore, we repeat the simulation now assuming a uni- 
form distribution of the radial velocity of RXJl 856.5-3754 in the 
range of Vr — —250 . . . + 250 km/s (see above discussion on the 
bow shock; note that we do not give priority tovr — O km/s). The 
results are summarised in Table 1 (b). Indeed, the parallax turns out 
to be larger (distance smaller) and the radial velocity tends to be 
close to zero. The derived values [Table 1 (b)] for the proper mo- 
tion and radial velocity imply an inclination of i = 88 ± 6 dcg, 
consistent with the observation of the bow shock (van Kerkwijk & 
Kulkami 2001). The displacement between RX J1856.5-3754 and 
the US centre would be 8.8 pc as inferred from the theoretically 
expected curve. The dmin distribution and the distribution of the 
corresponding flight times t are shown in Fig. 1, bottom panel. 
Assuming contemporary star formation, the mass of the progeni- 
tor star of RX J1856.5-3754 is given as the age of US (5 ± 2 Myr, 
e.g. Diehl et al. 2010; Preibisch et al. 2002, 2001; Preibisch & Zin- 
necker 1999; Blaauw 1978; de Geus et al. 1989) minus the time 
since the potential SN (~ 0.5 Myr). Using evolutionary models 
from Tinsley (1980), Maeder & Meynet (1989) and Kodama (1997) 
we estimate it to be 45 ± 3Mq, 43tg^ M© and 37+9^ Mq, re- 
spectively, corresponding to a main-sequence spectral type of 06, 
05 to 07 and 05 to 07 for the respective model (Schmidt-Kaler 
1982). Theprogenitor star of RXJl 856. 5-3754 should have an ear- 
lier spectral type than the earliest present member which has spec- 
tral type BO (HIP 81266, Madsen et al. 2002). Although NSs are be- 
lieved to form from progenitors with masses of < 20 — 25 Mq (e.g. 
Heger et al. 2003) it is known that in binary systems also more mas- 
sive stars, > 40 Mq, may produce NSs (e.g. Belczynski & Taam 
2008; Fryer et al. 2002; van den Heuvel & Habets 1984). Then, it 
should have ejected a runaway star, which may be known or un- 
known, yet. 



Table 1. Predicted current parameters of RXJl 856.5-3754 if it was born 
in US using as input for the radial velocity (a) a probability distribution 
according to Hobbs et al. (2005) and (b) a uniform distribution in the range 
of Vr = —250 . . . + 250 km/s; Vsp denotes the space velocity. The time 
and position of the SN is given as well. For the deduction of the radial 
velocity Vr, parallax vr, proper motion (fi^,fig), the distance of the SN 
from the Sun cIq (at the time of the SN), the SN equatorial coordinates 
(a, S, J2000.0, as seen from the Earth at present) and the time in the past 
at which the SN occurred (t, the kinematic age of the NS), see additional 
online supporting information. The distance of the SN to the US centre was 
obtained by adapting the theoretically expected curve (equation 1 in TIO) 



the slope of the dmin distribution. 
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Predicted SN position 
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Figure 1. Left: Distribution of minimum separations dmin between 
RXJ1856. 5-3754 and the US centre. Drawn as well is the theoretically ex- 
pected curve. Right: Distribution of corresponding flight times t in the past 
since the SN. 

Top panel: Distributions for Vr distribution derived from the Hobbs et al. 
(2005) distribution of pulsar space velocities. Solid curve in the left panel: 
fi = 7.9, (T = 2.0. Bottom panel: Distributions for Vr = —250 . . . + 
250 km/s. Solid curve in the left panel: fi = 8.8, a = 1.8. 
This indicates that the closest approach between RXJl 856.5-3754 and the 
US centre was 8 — 9 pc almost half a million years ago. 



3.2 Searching for a former companion of RX J1856.5-3754 

From section 3.1, US should be the parent association of 
RXJl 856.5-3754. Now, our attempt is to find a runaway star that 
might have been the former companion of the NS's progenitor 
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With US hosting the birth place of RX J 1856.5-3754, we know that 
the radial velocity of the NS is very small (Vr — 6^20 km/s). For 
the following analysis we adopt Ur = ± 50 km/s. 
Similar to section 3.1, we calculate the trajectories of RX J1856.5- 
3754 and any young runaway star with full 3D kinematics from 
Til, 1705 stars in total. 

First, we perform 10^ runs to find runaway stars which might have 
crossed the path of RX J 1856.5-3754 in the past. We find 16 stars 
for which the smallest separation to the NS found was less than 
10 pc (see additional online supporting information). Given their 
ages and kinematics, only three of them, HIP 74717, HIP 76028 
and HIP 78681, could have originated from the US association 
(0.5l2;iMyr, 1.8lE!:^Myr and 0.ll5!;i Myr ago, respectively) 
and only for them close encounters with RX J 1856.5-3754 are 
found to be situated inside US. After further three million runs 
for these three stars, we exclude HIP 76028 since the smallest 
separation found was larger than one parsec, see additional online 
supporting information for details on the selection criteria; not 
surprisingly since we find that the star probably originated from or 
at least passed US ^ 1 Myr before RX J 1856.5-3754 was born; if 
it is a BSS runaway star it may have formed in another SN in the 
Scorpius-Centaurus associations (14 to 20 SNe already exploded, 
Fuchs et al. 2009). 

For HIP 74717 the smallest dmin found was 0.06 pc. 
Adapting the theoretical distribution for 3D separations to the 
slope of the dmin distribution (restricting to 7558 runs where 
both, the runaway star and RX J 1856.5-3754 were within the US 
boundaries, Rus = 15 pc) implies a fly-by of the runaway and 
RX J 1 856.5-3754 with a displacement to each other of 7.3 ± 2.8 pc 
0.45loo2Myr in the past. Moreover, HIP 74717 is a single- 
lined spectroscopic binary with a rotational velocity of 38 km/s 
(Cutispoto et al. 2002), hence probably not a BSS runaway star. 
In the case of HIP 78681, the smallest dmin value found was 
0.04 pc. In only 329 Monte Carlo runs the distance of both, the 
runaway star and RXJ1856. 5-3754, to the US centre was less than 
15 pc, the nominal radius of the association. These runs yield dmin 
values (i.e. separations between the NS and the runaway star) from 
~ to ~ 25 pc. To improve statistics to be able to compare the 
dmin histogram with the theoretically expected curve, we include 
runs for which both stars were within 20 pc from the US centre, 
i.e. in a region at the US border The number of runs rises to 
49277. Then, the number of runs yielding large separations dmin 
is large; however, the peak at very small dmin values becomes 
well visible. For the first bins of this histogram, the theoretical 
curve suggests that both objects could have been very close to 
each other (positional difference ± 1.1 pc) 0.531qo5 Myr in the 
past. In the Simbad® database, HIP 78681 is listed with spectral 
type G7II whereas Pourbaix et al. (2004) list it as barium star 
(G8IIIBal). Moreover, HIP 78681 is a single-lined spectroscopic 
binary. Hence, the system is probably old with the companion 
of HIP 78681 being a white dwarf as widely excepted for SBl 
barium stars (Bohm-Vitense et al. 1984). However, if HIP 78681 is 
a young BSS runaway star, it might be the former companion to 
RX J 1856.5-3754 (although it is difficult to judge since relaxation 
of the US boundary was necessary). 

Thus, we do not find a convincing runaway star in the sample 
of runaway stars with full kinematics from Til to be a suitable 
former companion candidate for RXJl 856.5-3754. 



http://simbad.u-strasbg.fr/simbad/ 



Note, that the classical runaway star HIP 81377 (=(" Ophiuchi) is 
again not found to be the former companion to RXJ1856. 5-3754 
as it was suggested by Walter et al. (2000) but excluded by 
Hoogerwerf et al. (2001) and TIO (we found a smallest separation 
of 18.8 pc to the US centre after three million Monte Carlo runs). 

Since there are also 842 runaways in the runaway star cata- 
logue (Til) without radial velocities, we examined whether one 
of those could have been close to RXJ1856. 5-3754 in the past. 
We varied the radial velocity for those stars randomly within 
±500 km/s (the largest radial velocity values among all catalogue 
stars are ~ ±400 krn/s). After 10'' Monte Carlo runs 24 stars were 
found for which close encounters (^ 10 pc, see additional online 
supporting information) might have been possible in the past five 
million years. Nine of them showed dmin values smaller than 1 pc 
after 10' runs, hence were chosen for a more detailed investigation 
with 10® runs (for justification of the chosen limits see additional 
online supporting information). 



The distribution of the peculiar spatial velocities Vs 



of 



the population of young runaway stars is well represented by 
a Maxwellian distribution with a velocity dispersion of a = 
24.4 km/s (Til) and a maximum of 34.5 km/s. Extraordinary 
high velocities are unlikely although possible for individual cases. 
Only three stars - HIP 63803, HIP 70438 and HIP 74219 - 
would need plausible spatial velocities for close encounters with 
RXJl 856.5-3754; for the other stars, the necessary spatial veloc- 
ity deviates from the distribution maximum by more than 6a. In 
Table 2 the properties of the potential close encounter between 
RXJl 856.5-3754 and each of the three runaway stars are given. 
Considering the large encounter separation in the case of HIP 
70438, it is unlikely that HIP 70438 could have been at the same 
place as RX J1856.5-3754 in the past. 

In the case of HIP 63803 the position of the close encounter would 
have been far outside US. The spectral type of HIP 63803 is 
K2 (Simbad) and its position in the Hertzsprung-Russel diagram 
(HRD) suggests it to be a giant with an age of ~ 50 — 90 Myr 
(Til). Provided that a potential companion to that star formed at 
the same time and comparing to the proposed encounter time of 
^1.5 Myr in the past, the progenitor of RXJl 856.5-3754 should 
have been at least as old as 48 Myr when it exploded in a SN. Stars 
that end their lives in SNe do not live that long. Hence, HIP 63803 
cannot be the former companion to RXJl 856.5-3754 and might not 
be a BSS runaway star. 

In Tl 1 the F5V star HIP 74219 is fisted with a mass of 1.2±0. 1 Mq 
and an age of 39 ± 19 Myr that were determined from the median 
value of different evolutionary models not taking into account the 
error on the luminosity (mainly due to the parallax uncertainty) of 
the star. In the HRD the star lies just below the ZAMS of the evolu- 
tionary models. For that reason, it was treated as ZAMS star in Tl 1 . 
For higher luminosities according to its luminosity error the mod- 
els predict ages up to a few Gyr. Hence, it seems unlikely that HIP 
74219 was ejected from the young US association (~ 5 Myr), but 
we cannot fully exclude it as the former companion of RXJl 856.5- 
3754. 

So, we also do not find a convincing runaway star in the sample of 
runaway stars without radial velocity measurements from Tl 1 to be 
a suitable former companion candidate for RX J1856.5-3754. 



3.3 Identifying the parent association of RX J0720.4-3125 

Potential parent associations of RX J0720.4-3125 are investigated 
in the same way as for RXJ1856. 5-3754 (section 3.1). We repeat 
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Table 2. Predicted properties of potential close encounters between 
RXJ1856. 5-3754 and the three runaway stars without known radial veloc- 
ity measurements. The first column gives the Hipparcos identifier. Columns 
2 to 5 give the distance of the encounter to the centre of US (dus)> the 
separation d inferred from the theoretical curve (p it A), the time t of the 
encounter in the past as well as the peculiar spatial velocity Vsp,run of the 
runaway star needed. 



HIP 



c^US 
[pc] 



[pc] 



[Myr] 



[km/s] 



63803 240 .. . 270 
70438 50 ... 70 
74219 0...35 



0.0 ± 12.0 
6.7 ±3.4 
0.0 ±4.0 



,+0.15 -,70+95 

-0.10 -^''^-52 

Og.,+0.10 10+158 

-+0.08 cc+22 

-0.04 ""-IS 



1.42^ 



0.55^ 



the investigations already carried out in TIO here again as a new 
parallax measurement was done by Eisenbeiss (2011), yielding a 
distance of 280lgg" pc. Compared to the old value of 360lg™ pc 
(Kaplan et al. 2007), this new distance is in much better agreement 
with estimations derived from the spectrum and hydrogen column 
density uh (Posselt et al. 2007) giving 250 ± 25 pc although it is 
consistent within the error bars with the older value. 
In contrast to RX J1856.5-3754, the radial velocity of RX J0720.4- 
3125 is unconstrained (no bow shock detected) and is therefore 
derived from the probability distribution of pulsar space veloci- 
ties by Hobbs et al. (2005). We adopt the following parameters for 
RXJ0720.4-3125: 

Q = 07''20™24!961, 5 = -31°25'50"21 

(Kaplan et al. 2003), 
TT = 3.6 ± 1.6 mas (Eisenbeiss 2011), (2) 

/Iq = —92.8 ± 1.4 mas/yr (Eisenbeiss 2011), 
/if = 55.3 ± 1.7 mas/yr (Eisenbeiss 2011). 



For RXJ0720.4-3125, 18 associations/clusters are found for 
which separations dmin are consistent with the associations/cluster 
boundaries after 10** Monte Carlo runs. We then adapt equation 1 in 
TIO to the first bins of each dmin distribution to obtain the distance 
d of the SN to the association centre. Comparing the radii of each 
association with this putative separation d of the SN from the asso- 
ciation centre, nine associations/clusters are found to be potential 
birth places of RXJ0720.4-3125, i.e. d is consistent with the asso- 
ciation boundaries within its standard deviation: TWA, Tuc-Hor, /? 
Pic-Cap, the HD 141569 group, AB Dor, Collinder 140 (Col 140), 
Tr 10 and the Carina (CarA) and Argus Associations. In Table 3 
we give the position of the SN and the properties RX J0720.4-3 125 
would currently have if it was born in the respective association. 
Note that for the HD 141569 group, AB Dor and Tr 10 the proposed 
SN position is near the edge of the respective association. The HD 
141569 group may be excluded from the list of potential parent as- 
sociations since it contains only three stellar systems (five stars in 
total, HD 141569 is a triple system itself, Weinberger et al. 2001). 
Here, we retain it in the list of potential birth places (Table 3) for 
completeness. 

AB Dor currently contains young dwarfs with spectral types rang- 
ing from mid-F to early M (earliest members: HIP 18859 and HIP 
19183 with F5, Zuckerman et al. 2004; Lopez-Santiago et al. 2006) 
and masses from 0.7 to 1.3 Mq [derived from their positions in the 
HRD using evolutionary tracks (see Til for references)]. Hence, 
it seems unlikely that a SN occurred in ABDor although a 10 M© 
star might be expected to have formed in AB Dor from the compar- 
ison of its mass function with the initial mass function (Kroupa & 



Weidner 2005). 

Tr 10 is listed as sparse open cluster with a diameter of ~ 4pc 
(29' at 424 pc) in the Open Cluster Catalogue (Dias et al. 2002). 
However, de Zeeuw et al. (1999) and Brown et al. (1999) found 
it to be actually an association with a diameter of ~ 50 pc (8° at 
366 pc) rather than a cluster. The 23 members found by de Zeeuw 
et al. (1999) are all B3 to AO type stars. Since early B type stars are 
present in Tr 10, it is plausible that the association already experi- 
enced a SN. 

Close encounters between RX J0720.4-3125 and Col 140 occur at 
0.04 ± 0.01 Myr in the past. In this case, the SNR should still be 
visible; however there is no known SNR in this area (Green 2009; 
A. Poghosyan, priv. comm.). 

Hence, except the HD 141569 group and the open cluster Col 140, 
the associations listed in Table 3 are still probable birth associations 
of RX 10720.4-3 125 although AB Dor seems unlikely. 

3.4 Searching for a former companion of RX J0720.4-3125 

After lO"* Monte Carlo runs 53 runaway stars with full kinematics 
show a smallest separation dmin to RX J0720.4-3125 of less than 
10 pc. For them, another three million runs are performed (see 
additional online supporting information). We select those stars 
for which separations smaller than one parsec are found after three 
million runs (see additional supporting information), 29 in total. 
To further reduce the number of potential former companion 
candidates, we investigate whether the potential close encoun- 
ters occurred inside one of the possible birth associations of 
RX J0720.4-3125 listed in Table 3. Possible close encounters 
between RX J0720.4-3125 and 16 stars are found to be consistent 
with the radius of one or more of the nine associations/clusters. 
Of them, we can exclude HIP 36160 (spectral type G5V) whose 
chromospheric age is 3.27 Gyr (Rocha-Pinto et al. 2004), the 
cataclysmic variable HIP 40430 (spectral type 09nne) and the 
barium star HIP 47267 (Bergeat & Knapik 1997; spectral type 
G8II). The latter has been included in the sample of Til owing 
to its luminosity class. We also exclude HIP 63972 (spectral type 
KOII-III) that has been classified as old disk star (Eggen 1993) and 
entered the sample of Til due to its luminosity class. The ages 
for HIP 67655 (spectral type G5V), HIP 68783 (G5), HIP 74689 
(spectral type A4V) and HIP 79377 (spectral type GI/G2V) are 
very uncertain since these stars currently lie on or close to the 
main-sequence. From its position in the HRD we estimate the age 
of HIP 74689 to be roughly 50 to 300 Myr, too old for an origin 
in P Pic-Cap (age 8 — 34 Myr) which is the association found 
to have hosted the potential encounter. The other three stars lie 
slightly below the model zero-age main sequence (ZAMS), thus 
were treated as ZAMS stars in Til. For HIP 79377 Holmberg 
et al. (2009) give an age of 200 Myr. With a later spectral type, 
HIP 67655 and HIP 68783 are probably older. For those reasons, 
we do not consider HIP 67655, HIP 68783, HIP 74680 and HIP 
79377 as former companion candidates for RXI0720.4-3125. 
At this stage, we do not exclude binary (and multiple) runaway 
systems as potential former companion candidates to the NS 
progenitor. The reason is that it is not excluded that a former binary 
(in rare cases also multiple) companion could have survived the 
SN explosion. 

For the remaining eight stars - HIP 40326, HIP 43158, HIP 57269, 
HIP 59803, HIP 76304, HIP 78078, HIP 78106 and HIP 84794 
- we explore the distribution of dmin- We adapt the theoretically 
expected curve (equation 1 in TIO) to the first part of the dmin 
distribution. 
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Table 3. Potential parent associations of RXJ0720.4-3125. 

In the second and third Columns we give the predicted encounter separation d (inferred from the theoretical curve) and encounter time t (determined from 
those runs which yield separations consistent with d). Columns 4 to 8 give the predicted present NS parameters (radial velocity Vr, proper motion /i* and fig, 
parallax n) and space velocity Vsp for each case and Columns 9 to 1 1 indicate the distance to the Sun d^ (at the time of the SN) and equatorial coordinates 
(J2000.0, as seen from Earth at present) of the potential SN. Error bars denote 68 per cent confidence (cf. Appendix B, TIO). 
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Table 4. Potential former companion stars of RXJ0720.4-3125 and the respective parent association/cluster. For details on the determination of the predicted 
NS parameters given in Columns 4-8 (radial velocity Vr, proper motion ^* and fig, space velocity Vsp, parallax w) see additional online supporting informa- 
tion. Columns 9-11 indicate the distance to the Sun i^q and equatorial coordinates (J2000.0, as seen from the Earth at present) of the potential SN (at the time 
of the SN). The last Column gives an estimation of the mass of the progenitor star derived from the progenitor's life-time (age of the runaway star, see text, 
minus the time since the potential SN) using evolutionary models from Tinsley (1980), Maeder & Meynet (1989) and Kodama (1997). 

For HIP 43158 the results using a Vr distribution derived from the one of pulsar spatial velocities (Hobbs et al. 2005) (superscript *) as well as a uniform 
distribution in the range Vr = —300 . . . + 300 Icm/s (superscript #) are shown (see text). 
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Three stars appear to have possibly been at the same place as 
RXJ0720.4-3125 (i.e. fj, = 0): HIP 43158, HIP 57269 and HIP 
76304. In Table 4 we give the time and position of the SN and the 
properties RX J0720.4-3125 would currently have if it was born in 
the respective association with the runaway star being the former 
companion. In the last column of that table an estimation of the 
mass of the progenitor star is given derived from the age of the 
runaway star (see below) minus the time r since the potential SN 
using evolutionary models from Tinsley (1980), Maeder & Meynet 
(1989) and Kodama (1997). 

In the case of HIP 43 158 we find the radial velocity of RX J0720.4- 
3125 to be rather small. For that reason, we repeat the calculations 
adopting a uniform radial velocity distribution in the range of 
-300 . . . + 3001^m/s (cf. RXJ1856.5-3754 and US, section 3.1). 
For HIP 43158, the results for both radial velocity distributions are 
given in Table 4. 

HIP 43158 is a single BO giant star with a peculiar space ve- 
locity of ~ 53 km/s and an age of ~ 15 Myr (Til) which is 
consistent with an origin in Tr 10 (age 15 — 35 Myr, see TlO 
and references therein). It has previously been suggested that 
HIP 43158 originates from the Vela region (Schilbach & Roser 
2008). The position of the proposed SN (from uniform Vr, case 
#: I = 258.9 . . . 259.9 deg, b = 1.8 .. . 2.5 deg) (that occurred 
« 0.9 Myr ago) lies in a region of enhanced ^'^ Al emission (Diehl 



et al. 2010, their figure 2, ^''AI lifetime ^ 1 Myr). Since ^''AI 
is produced by Wolf-Rayet winds as well as SN explosions (for 
a review see Prantzos & Diehl 1996), we expect to find ^''Al 
emission at the position of a recent SN. 

HIP 57269 (spectral type KIV) is a double star in a triple system 
(Konig et al. 2003) with a peculiar space velocity of ~ 21 km/s 
and an age of ~ 20 Myr (Til) that is still in its pre-main sequence 
phase (Konig et al. 2003). Its age is consistent with that of Tuc-Hor 
(10 — 40 Myr, see TIO and references therein), the potential host 
association of the SN. Comparing the position of the proposed SN 
(l = 281.3 . . . 291.1 deg, b = 15.3 . . . 20.6 deg) (that occurred 
^ 0.3 Myr ago) with the map of ^''AI 7-rays (Diehl et al. 2010), 
there seems to be only little '^''Al emission in this area. Also, it 
might be doubtful whether a triple system could have survived a 
SN in a multiple system. 

The suspected spectroscopic binary HIP 76304 (Frankowski et al. 
2007; spectral type G2V) is an X-ray source and is listed as T 
Tauri star in the catalogue of T Tauri stars in Scorpius-Centaurus 
(Sco-Cen) by Kohler et al. (2000) but with unknown classification. 
It has a peculiar space velocity of ~ 45 krri/s. Holmberg et al. 
(2009) give an age of 3.5 to 5.9 Gyr which would be far too old to 
be associated with Sco-Cen or another young nearby association. 
Since it is an X-ray source, a pre-main sequence star with an age of 
~ 9 Myr (Til) seems more plausible. Then its age would be con- 
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Figure 2. (B — V)q versus My diagram of Tr 10 (dots represent mem- 
ber stars according to de Zeeuw et al. 1999, the star marJis HIP 43158). 
The solid and dashed lines represent the 15 Myr and 35 Myr isochrones 
from Marigo et al. (2008) (for solar metallicity; http://stev.oapd.inaf.it/cgi- 
bin/cmd). {B — V)q and My of the stars are derived from their spectral 
types according to Schmidt-Kaler (1982). 



sistent with that of /? Pic-Cap (8 — 34 Myr, see TIO and references 
therein), the potential host association of the SN. At the position of 
the supposed SN (l = 336.1 . . . 341.4 deg, b = 31.5 . . . 33.2 deg) 
(that occurred ~ 0.6 Myr ago), an enhancement of '^^Al emission 
(Diehl et al. 2010) seems to be present. 

The rotational velocities («sini) of these three stars are rather 
small; vsini = 96 ± 15km/s, 20km/s and 5km/s for HIP 
43158, HIP 57269 and HIP 76304, respectively (Penny 1996; 
Randich et al. 1993; Holmberg et al. 2009). If one of them 
originated from a supernova in a multiple system, this may indicate 
a small inclination angle i. 

Considering that in the case of HIP 57269 only little ■^'' Al emission 
is visible although the potential SN should have occurred very 
close to the Sun (~ 27 pc) and only ~ 0.3 Myr ago, this scenario 
seems less likely than for the other two candidate stars. Looking 
at the proposed parallax RX J0720.4-3125 would currently have 
in each case, also HIP 76304 is a less good candidate since the 
parallax would be much larger (5.6 — 6.8 mas) than the measured 
value (3.6 ± 1.6mas, Eisenbeiss 2011). 

We suggest that RXJ0720.4-3125 may have been born in a SN 
0.85 ± 0.15 Myr ago as a former member of Tr 10 with HIP 
43158 being the possible former companion. It has been suggested 
that BSS runaway stars should be blue stragglers due to mass 
transfer during binary evolution, i.e. they appear younger, hence 
bluer, than their parent association (see also Hoogerwerf et al. 
2001 for other examples). In Fig. 2, we show the positions of HIP 
43158 and Tr 10 member stars from de Zeeuw et al. (1999) along 
with isochrones for 15 Myr and 35 Myr taken from Marigo et al. 
(2008). If Tr 10 is only as young as 15 Myr, HIP 43158 is not 
a blue straggler; however, Tr 10 can be as old as 35 Myr (TIO 
and references therein). In this case HIP 43158 would be a blue 
straggler. 

The mass of the progenitor star would have been 13 — 14 M© 
(for an age of Tr 10 of 15 Myr) to 7 — 9 Mq (for an age of Tr 
10 of 35 Myr) corresponding to a spectral type of Bl to B2/3 on 
the main-sequence. This is consistent with the progenitor star of 
RXJ0720.4-3125 having an earlier spectral type than the earliest 
current member of Tr 10 (four B3 stars, de Zeeuw et al. 1999). 
Tr 10 has been previously suggested to host the birth place of 



RXJ0720.4-3125 by Motch et al. (2003) who considered the 
general direction of the NS's motion and Kaplan et al. (2007) who 
investigated the probability of close approaches of the NS to any of 
the OB associations given in de Zeeuw et al. (1999). They varied 
the parallax within 2.8 ± 0.9 mas and the radial velocity in the 
range Vr = ±0.935i't (0.935 corresponds to la in Vr for random 
orientation, Vt is the transverse velocity) and found a separation 
between the NS and the centre of Tr 10 of 17 pc 0.7 Myr ago for 
Vr = —20 . . . + 50 km/s, not inconsistent with our more complete 
calculations. 

Since we cannot pre-constrain the radial velocity of 
RX J0720.4-3125, it is difficult to examine those 842 runaway star 
candidates in Tl 1 without radial velocity measurements due to the 
large uncertainties (assuming Vr = —500 . . . 500 km/s). For al- 
ready 304 of the 842 stars it is possible to find a past position as 
close as 10 pc to RX J0720.4-3 125 after lO'' Monte Carlo runs. Af- 
ter further three million runs, we found for 1 13 of those 304 stars 
a smallest separation to RXJ0720.4-3125 that was less than one 
parsec (for justification of the limits see additional online support- 
ing information). Excluding stars that would need a peculiar space 
velocity larger than 180 km/s (that is ~ 6a above the maximum 
of the distribution of runaway star velocities, Tl 1), 34 candidates 
are left. For all of them we can find a set of input parameters (vr, 
H^, fr) for that we can find close encounters (dmin ^ 10 pc) with 
RX J0720.4-3125 within one or more of the nine potential birth 
places listed in Table 3. Here, constraining the runaway star radial 
velocity through observations is necessary. 

Most of the 34 stars are late-type stars that are situated in the HRD 
on or near the ZAMS. For that reason, we doubt that they are in- 
deed sufficiently young runaway stars for RX J0720.4-3125 (< few 
Myr).^ 



4 SUMMARY AND CONCLUSIONS 

We analysed the origin of two members of the M7, RX J 1856. 5- 
3754 and RX J0720.4-3125 using most recent parallax measure- 
ments. 

Under the assumption that most NSs are born in and ejected from 
their parent association or cluster (see section 1), we confirm that 
RX J 1856.5-3754 most probably originated from the US associa- 
tion as suggested previously (Walter & Lattimer 2002, TIO). We 
find that the radial velocity of RX J1856.5-3754 is GtUkm/s 
implying an inclination angle to the line-of-sight of 88 ± 6 deg that 
is consistent with the inclination of the bow shock that RX J 1 856.5- 
3754 creates in the ISM owing to its motion (van Kerkwijk & 
Kulkami 2001). This consistency implies that the current distance 
of RX Jl 856.5-3754 is « 120 pc as obtained by Walter & Lattimer 
(2002) and Walter et al. (2010) rather than 160 - 180 pc (Kaplan 
2003; van Kerkwijk & Kaplan 2007). This smaller distance is also 
in good agreement with RXJl 856.5-3754 being in front of the 
Corona Australis star forming region (~ 130 pc, Neuhauser & 
Forbrich 2008); if it was behind, optical detection would probably 
have failed. Moreover, this is in accordance with the absence of OI 
from the ISM in the Chandra HRC/LETG spectra (Burwitz et al. 
2001, 2003). The space velocity of RXJl 856.5-3754 would then 
be ~ 195 km/s. The derived kinematic age of RXJ1856. 5-3754 

^ In Tl 1, they were treated as ZAMS stars and contribute to the 20 per cent 
contamination that are no true runaway stars but fast-moving members of 
the low-velocity group. 
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of 0.46 ± 0.05 Myr agrees well with theoretical cooling models 
(cf. TIG, Figure 12 therein). 

The mass of the progenitor star would have been ~ 40 — 60 M0. 
Considering that the NS may have originated from a binary 
system that got disrupted in an asymmetric SN explosion, we 
searched for a possible former companion that should then be a 
runaway star. No former companion candidate was identified. A 
reason could be that, if there was a former companion, either, it 
may not be included in the runaway star catalogue of Til either 
because it was no Hipparcos source or its velocity vector is not 
significantly different from those of its neighbouring stars (note, 
that Til also investigated the direction of motion of young stars 
compared to neighbouring stars, regardless whether the absolute 
velocity is large or small); or a former companion might also 
have experienced a SN after ejection (then, the mass ratio of 
the binary would have been close to unity). Alternatively, the 
progenitor star of RX J 1856.5-3754 could have been a single 
star and the NS gained its high velocity in an asymmetric SN 
explosion. However, single stars of masses > 20 — 25 M© 
are generally believed to produce black holes rather than NSs. 
Maiz-Apellaniz (2001), Berghofer & Breitschwerdt (2002) and 
Fuchs et al. (2009) suggested that about 20 nearby (in Sco-Cen) 
SNe created and re-heated the Local Bubble. RXJl 856.5-3754 
was probably formed in one of them about half a million years 
ago in US. For that reason and since the small radial velocity 
of RXJ1856. 5-3754 that we found is strongly supported by the 
bow shock measurements of van Kerkwijk & Kulkarni (2001), 
we argue that RXJl 856.5-3754 is not a remnant of a runaway 
star that had already left its parent cluster when it experienced a SN. 



temperature of RX J0720.4-3 125 is fcT ^ 85 - 95 cV (Hohle et al. 
2009) while for RXJl 856.5-3754 it is kT ^ 60 eV (Burwitz et al. 
2003). Hence, RXJl 856.5-3754 is cooler although it is probably 
younger. This might be due to different birth parameters or differ- 
ent cooling. Moreover, effective temperatures may be influenced 
by hot spots, hence do not reflect the real surface temperature. For 
RX J0720.4-3125, this is indeed the case (Haberl et al. 2006; Hohle 
et al. 2009) and it appears hotter. 
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For RXJ0720.4-3125 there is no unique result on the parent 
association due to the large parallax uncertainties (as long as run- 
away stars are not concerned). It could have been formed either in 
one of the young local associations (Fernandez et al. 2008) with 
TWA being the best candidate birth place (see also TIO) or in a 
more distant association/cluster such as Tr 10 or Col 140. In the first 
case, the distance of the SN to the Sun is very small (~ 30 — 50 pc) 
and it is expected to find SN-produced radionuclides on Earth (Ellis 
et al. 1996). A small but insignificant signal of ''"Fe was found for 
0-1 Myr in the past (Knie et al. 2004). 

The identification of another indicator is thus needed. We find 
three runaway stars for which not only a positional coincidence 
with RX J0720.4-3125 is possible but also an enhancement of '^^Al 
emission was measured at the location of the potential SN (Diehl 
et al. 2010). Regarding the current distance of RX J0720.4-3125 
derived for each case, we suggest that RX J0720.4-3125 may have 
been born 0.85 ± 0.15 Myr ago as a former member of Tr 10 
with HIP 43158 being the possible former companion. The mass 
of the progenitor star would have been ^ 8 — 14 M©. The cur- 
rent distance of RXJ0720.4-3125 would then be ^ 285 pc, in 
good agreement with estimates from nn measurements (250 ± 
25 pc, Posselt et al. 2007) and the most recent HST parallax (dis- 
tance 280^85" pc, Eisenbeiss 201 1). The present radial velocity of 
RX J0720.4-3125 would be -76^17 km/s leading to a space ve- 
locity of ~ 160 km/s. 

However, we emphasise that further observation of the former com- 
panion candidate is needed to confirm its status. There are also run- 
away stars without radial velocity measurements that are potential 
former companion candidates for RX J0720.4-3125. For them, a 
measurement of the radial velocity is crucial. 
RXJ0720.4-3125 is probably older than RXJl 856.5-3754, thus 
should be cooler than RXJ1856. 5-3754. However, the effective 
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